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Abstract 



Test modes on turbulent magnetized plasmas are studied taking 
into account the ion trapping that characterizes the E x B drift in the 
background turbulence. We show that trappyng provides the phys- 
ical mechanism for the formation of large scale potential structures 
(inverse cascade) observed in drift turbulence. Trapping combined 
with the motion of the potential with the diamagnetic velocity deter- 
mines ion flows in opposite directions, which reduce the growth rate 
and eventually damps the drift modes. It also determine transitory 
zonal flow modes in connection with compressibility effect due to the 
polarization drift in the background turbulence. 

Keywords: plasma turbulence, nonlinear processes, structure gener- 
ation. 



1 Introduction 

The evolution of turbulence in magnetically confined plasmas is a complex 
problem that is not yet completely understood besides the huge amount of 
work on this topic (pQ and the references there in). Low-frequency drift 
type turbulence, which has significant influence on the magnetic confinement 
of high temperature plasmas, is extensively studied especially in connection 
with fusion research (see e.g. [2], [3], jl]). Most of the studies that go beyond 
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the quasilinear stage are based on numerical simulations or on simplified 
models. They show a complex nonlinear evolution with generation of large 
scale structures, increase of order and appearance of zonal flow modes ([5], 
[6]) that leads to the nonlinear damping of turbulence. 

The stochastic particle advection that appears in turbulent plasmas due 
to the E x B (or electric) drift determined by the fluctuating potential can 
produce trajectory trapping or eddying around contour lines of the potential. 
Trapping has a strong effect on the statistics of particle trajectories. Ana- 
lytical methods adequate to the study of particle stochastic advection in the 
presence of trapping were developed only in the last decade ([7], [8] and the 
references therein). They permitted to understand that trapping determines 
strong departure from the characteristics of Gaussian advection processes. 
The aim of this paper is to contribute to the understanding of the effects 
of trajectory trapping on the evolution of drift turbulence. These are the 
first analytical results on this complex problem that are in agreement with 
numerical simulations. 

Test particle trajectories are strongly related to plasma turbulence. Plasma 
dynamics is basically described by the Vlasov-Maxwell system of equations, 
which represents the conservation laws for the electron and ion distribution 
functions along particle trajectories coupled with the constraints imposed by 
Maxwell equations. Analytical studies of plasma turbulence based on tra- 
jectories were initiated by Dupree [9], [TO] and developed especially in the 
years seventies PQ. These methods do not account for trajectory trapping 
and thus they apply to the quasilinear regime or to unmagnetized plasmas. 
A very important problem that has to be understood is the effect of the 
non-standard statistical characteristics of the test particle trajectories on the 
evolution of turbulence in magnetized plasmas. A Lagrangian approach is 
developed, which extends the Lagrangian methods of the type of [TO], [TT] . 
[T2] to the nonlinear regime characterized by trapping. 

We study linear modes on turbulent plasma with the statistical character- 
istics of the potential considered known. Drift turbulence in hot magnetized 
plasmas is considered. Analytical expressions are derived, which approxi- 
mate the growth rates and the frequencies of the test modes as functions of 
the characteristics of the background turbulence. They provide an image of 
turbulence evolution. 

We show that there is a sequence of processes, which appear at different 
stages of evolution as transitory effects and that the drift turbulence has 
an oscillatory (intermittent) evolution. A different perspective on important 
aspects of the physics of drift type turbulence in the strongly non-linear 
regime is deduced. The main role in these processes is shown to be played 
by ion trapping. 
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The paper is organized as follows. The dispersion relation for the drift 
modes on turbulent plasmas is deduced in Section 2. It is shown that the 
effects of the background turbulence are contained in a function of time that 
enters in the ion propagator. It is an average over the stochastic ion trajec- 
tories in the background turbulence. The statistical methods for evaluating 
this function are shortly presented in Section 3. The next three Sections 
present the effects of the turbulence on the test modes for quasilinear tur- 
bulence, weak nonlinear regime (when the fraction of trapped ions is small) 
and strong nonlinear regime (with strong trapping). The evolution of the 
drift turbulence is discussed at each stage. The summary of results and the 
conclusions are presented in Section 8. 

2 Test modes in turbulent plasmas 

Drift waves and instabilities are low-frequency modes generated in non- 
uniform magnetically confined plasmas. Depending on the particular con- 
ditions, there are several types of drift modes. Since the aim of this work 
is to understand the effects of trapping on the evolution of turbulence, we 
consider a simple confining geometry, the plane plasma slab, in which the 
magnetic field is straight and uniform. Plasma has low j3, which means that 
the perturbation of the magnetic field is negligible (electrostatic approxima- 
tion). 

The magnetic field is along z axis (B = Be z ) and plasma is non-uniform in 
one direction taken along x axis. For simplicity, the equilibrium temperatures 
are uniform and only the density no(x) is x-dependent. The characteristic 
length of density variation L n = n / \dno/dx\ is much larger than the wave 
length of the drift modes. 

We start from the basic description of this (universal) drift turbulence 
provided by the drift kinetic equation in the collisionless limit. Electron 
kinetic effects produce the dissipation mechanism to release the energy and, 
combined with finite Larmor of the ions, make drift waves unstable. The 
latter consist of the polarization drift velocity and of the modification of the 
electric drift velocity due to the gyro-average of the potential on the ion 
orbits. Both effects determine the decrease of mode frequency below the 
diamagnetic frequency, which makes the growth rate positive. Beside this, 
the polarization drift has a more complex influence determined by its nonzero 
divergence. We neglect here the modification of the potential and consider 
the effects of the polarization drift. The reason is that, as shown below, 
the background turbulence in the nonlinear regime eliminates the effects of 
the finite Larmor radius on the frequency. Thus, this approximation that 
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determines sensible simplifications of the calculations has negligible influence 
on the phenomenology of drift turbulence evolution in the nonlinear regime. 

The drift kinetic equations for the perturbation 5f a of the distribution 
function of the guiding centers / a (x, z, v z , t) = F a {y z ) n (x) + 5/ Q (x, z, v z , t) 
for electrons and ions is 

d t Sf a + u x • V5f a + v z d z 5f a - ^ (d z </>) d v Jf a = F a (v z ) d y d x n (1) 

m a 

where a represents the species (a = e, i), (ft is the potential, x = (x,y), 
F a (v z ) is the Maxwell distribution, v z is the velocity along the magnetic 
field and e a , m a are the charge and mass of the particles. The perpendicular 
velocity is ui = — V0 x e z /B. 

The perturbed densities 5n a are obtained by integrating the distribution 
function 5f a over the velocities v z and Vj_. Poisson equations, which is 
approximated by the quasi-neutrality condition Sn e = 5n l , closes the system. 

The instabilities are usually studied on quiescent plasmas by introducing 
in Eq. (jTJ a wave type potential <fi = 5<fi where 

S(j)(x, y, z, t) = <fi kuj exp {ik x x + ik y y + ik z z — iut) . (2) 

Linearizing the equations, the frequency u and the growth rate 7 are deter- 
mined from the quasineutrality condition, which is the dispersion relation for 
the drift waves. 

This model of modes developing on quiescent plasma is not realistic be- 
cause drift instabilities appear for a large range of wave numbers and produce 
a turbulent potential. Test mode models consider a turbulent plasma with 
given statistical characteristics of the background potential 0&(x, z,t) and a 
small perturbation 5(f>, (ft = <pb + 84>- The growth rates and the frequencies of 
the test modes are determined as functions of the statistical characteristics 
of the background potential 4>b by linearizing Eq. (pQ) around the background 
potential. 

The aim of the present study is to determine the effects of the turbulence 
on test modes and in particular the influence of trajectory trapping or eddy- 
ing produced by the electric drift. The dispersion relation of the drift waves 
in quiescent plasma (<pb = 0) is review in subsection 2.1 in order to have a 
comparison basis for the effects appearing in turbulent plasmas. The disper- 
sion relation for turbulent plasmas is determined in subsection 2.2, where we 
show that the background potential modifies the propagator of the modes 
through trajectory distribution and also by a compressibility effect produced 
by the polarization drift. 
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2.1 Drift modes in quiescent plasmas 

We begin with the electrons and show that their response is the same in quies- 
cent and turbulent plasmas because actually they do not "see" the turbulence 
due to the fast decorrelation produced by the motion along the magnetic field. 

Electrons are dominated by the parallel motion for the small frequency 
fluctuations with u <C k z VT e - The density of electrons can be written as 

/ e (x, z, v, t) = n Q {x)F« expA + h e (3) 

■L e 

where the first term is the adiabatic response obtained from the parallel terms 
in Eq. (pQ) and the non-adiabatic term h e is the solution of 

t h e - Yl^l . Vh e + vAh e _ f« (^0) d Sf a (4) 

B m a 
= n F e M (V* e d y + d t )^-, 

J- e 

where the diamagnetic velocity 

V - Te dn ° - c — 
en B dx\ s L n 

was introduced and the condition e<p/T e < 1 that holds for the drift waves 
was used. T e is the electron temperature, c s = ^jTjmi, p s = c s /VLi is the 
Larmor radius and f2j = eB/rrii is ion cyclotron frequency. 

Due to the large electron velocity v z , the perpendicular drift can be ne- 
glected in electron trajectories as well as the parallel acceleration which pro- 
duces velocity variations that are negligible compared to the thermal velocity. 
The trajectories are not influenced by the potential and consequently electron 
equation for drift waves can be approximated by the linear one 

d t h e + v z d z h e = n F e M (V* e d y + d t ) ^. (5) 

The solution is obtained by using the method of characteristics. The right 
hand side of Eq. (j^J) is integrated along the trajectory, which is approximated 
by z{r) = z — v z (t — r). For the potential (j2J) it is 

, e C0fcw ™ UJ — kyV* e . . . . . 

n = nob M — : exp uk.x + %k z z — tut) . (o) 

T e v z k z - cu 

The perturbation of the electron density is obtained by integrating f e over 
velocities. The integral aver v z , which is singular, is determined in the com- 
plex u plane by the pole of the singularity (the principal value is negligible) 
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[see, for example, [IB] , page 457]. One obtains 



5rf = n (x) 



7T oj — k v V, 



y * *e 



exp (ik.x + ik z z — iut) , (7) 



T e V V 2 |fc 2 |fT, 
which hold for both quiescent and turbulent plasmas. 

The ion response is influenced by finite Larmor radius effects: the po- 
larization drift determined by the time variation of the electric field and 
the modification of the potential by the gyro-average. These effects com- 
bined with the non-adiabatic response of the electrons destabilize the drift 
waves. In order to concentrate on the nonlinear effects produced by trajec- 
tory trapping, we make an additional simplification, which has the advantage 
of simplifying the analytical expressions: we neglect the change of the poten- 
tial by the gyro-average and take the same potential in the ion and electron 
equations. The reason is that the turbulence evolves in the nonlinear stage 
to correlation lengths much larger than the ion Larmor radius and the effect 
of potential gyro-average becomes negligible. As we are interested in this 
paper by effects that appear beyond the quasilinear stage, this simplification 
does not affect the results. 
The polarization drift 

:dtE± 



u, 



eB 2 



is much smaller than the electric drift (by a factor oj/Qi, where) and actually 
it has negligible effect on ion trajectories. The polarization drift is important 
in the ion equation due to its divergence 



u, 



eB 2 



d t A(/>, 



(9) 



which leads to the perturbation of the ion density. We consider perturbation 
with small parallel wave numbers u/k z ^> vn and the parallel motion can be 
neglected in Eq. ([1]) for the ions. The linearized equation with the potential 
© is 

eF i 

d t Sf = -no-zr-VndySt - n F l M V ■ u p (10) 



The solution is 



6f 



x,t, v) 



-no 



eF l M 



56 



kyV*e „ ojk i 



e 2 B 2 



IT 



(11) 



where k± = ^Jk 2 + k 2 . The propagator II 1 is defined by 



IT 



drexp (ik- (x(r) — x) (r) — iu(r — t)) 



(12) 



6 



with the integral taken along ion trajectories. In the linearized case II* = 
—1/oj. The perturbation of the ion density is obtained by integrating over 
velocities 

5n l (x, t, v) = n — [kyVte - uk±pl] . (13) 

J-e W 

The quasineutrality condition leads to the dispersion equation 

and to the well known solution for drift instability 

uj = -^5-, (15) 
l + k ± 

7 = jouJ (k y - u) , (16) 



where fc< = kip SJ u = coL n /c SJ 7 = ^L n /c s , 70 = ^/7r/2(c s /L n )/ \k z \ v Te . 

As seen in the above equations, the drift instability is determined by the 
nonadiabatic response of the electrons, which leads to 7 > if uj < k y V* e , 
condition ensured by the finite Larmor radius of the ions. This conditions is 
fulfilled by all values of k± and thus all the modes are unstable. The wave 
number domain of unstable modes is very large, and the maximum growth 
rate is for u = k y V* e /2, which corresponds to k±p s = 1. The growth rate is 
quadratic in k y and thus the most unstable modes have the largest values of 
k y compatible with the above condition. These are the characteristics of the 
linear (universal) drift instability on quiescent plasmas. 

The solution in the zero Larmor radius limit is u = k y V* e , 7 = 0, which 
represents the stable drift waves. For an arbitrary initial condition <f) Q , this 
solution is 

<f)(x, y, z, t) = <f) (x, y-Vtet, z). (17) 

Thus, the basic effect produced by the plasma when a potential appears 
is to displace it with the diamagnetic velocity. 



2.2 Dispersion relation in turbulent plasma 

We consider a turbulent plasma with given statistical characteristics of the 
stochastic potential 0&(x, t). The potential is taken as the zero order solu- 
tion (j!7p . obtained when the polarization drift is neglected. It consists of 
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the motion of the potential with the diamagnetic velocity. The modifica- 
tion of potential shape and amplitude appear due to polarization drift on a 
larger time scale of the order I/7. This approximation is confirmed by the 
numerical simulations (see for instance p3] where the Eulerian correlation 
of the potential obtained from numerical simulation of the trapped electron 
mode turbulence evidences the motion of the potential with the diamagnetic 
velocity). The test mode studies of turbulence are based on this time scale 
separation, which permits a sequential approach. Starting from a potential 
that is a zero order solution ( ITT)) it is possible to determine the frequency and 
the growth rate of test modes as function of the statistical characteristics of 
the potential. They provide information on the tendency in the evolution of 
the potential, which is used to determine the test mode properties later in 
the evolution, and so on. 

The main statistical characteristics of the background turbulence are the 
amplitude (3 of the potential fluctuations, their correlation lengths A^, \ y 
and correlation time r c . These parameters appear in the Eulerian correlation 
(EC) of the potential defined by 

E(pc, t) = (0 b (x', 0&(x' + x, t' + t)> , (18) 

where () is the statistical average or the space average. This function is the 
Fourier transform of the spectrum. The amplitude of the stochastic electric 
drift is V = -JV% + V y 2 , where V x = (3/B\ y , V y = (3/B\ x . These parameters 
define the time of flight (or the eddying time) r/j = X x /V x = \ x \ y B/(3, which 
is the characteristic time for trajectory trapping. 

The electron response to a perturbation with 5<fi of the background po- 
tential (fib is given by Eq. (jTJ), as shown in Section 2.1. 

In order to determine ion response, we define the operator of derivation 
along ion trajectories (where the parallel motion is neglected) 

, sa( _V^xb y (lg) 
The equation for the distribution function is 

af + f v • up = 0. (20) 

The distribution 

,i / \h A , e0 6 (x-V* e t)\ 

f = n Q {x)F l M (1 + K — >-\ (21) 

represents the approximate equilibrium because (7/o = and the term 
V ■ Up 1 is of the order V/\fliT c . The divergence of the polarization drift 
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has the dimension of t~ l and it introduces a characteristic time, which is of 
the order r p = r c {fljjfi) thus much larger than the correlation time of the 
potential. It means that for time much smaller than this the remaining term 
is negligible and the ion distribution function can be approximated by ( )2~Tj) . 
Perturbing the potential with 5(f), the operator is perturbed by 

5 Ct = (l/B) (-V5(j) xe 2 )-V 

and a change of the distribution function appears f l = fo + h. The linearized 
equation in this perturbation is 

(jh + Sdfo + hV -u p + f*V-5u p = 0. (22) 

Since the background potential remains small compared with the kinetic 
energy e0&/T e << 1, the equilibrium distribution function in the second and 
forth term can be approximated by n (x)Fl I 

&h + hV ■ u p = -in {x)Fi I ^- (k y V* e - up 2 s kl) . (23) 

This equations shows that the right side term is not influenced by the tur- 
bulence (it is the same as for quiescent plasmas (|T0|) ). The effects of the 
background potential appear in the trajectories (in the operator of derivation 
along trajectories) and in the second term which accounts for the divergence 
of the polarization drift produced by the turbulence. 
The formal solution is 

/i(x, v, t) = - no (x)F* M ^ (k y V* e - up 2 s kl) IT (24) 

J- e 

where the propagator is 

IP = % I <ir exp [ik- (x(r) — x) — icu (r — t)] exp — / V • u p (x(r)) 

J —OO L J T 



and the integrals are along ion trajectories obtained from 

dx(r) V0(x - V* e t) x e z 
dr ~ B ' 



(25) 
(26) 



with the condition at r = t, x(t) = x. 

The ion response (the propagator) is averaged over the stochastic trajec- 
tories 



if = i f dr M(r; t) exp [-iu (r - t)} 

J — oo 



(27) 
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where M is the average 



M(r; t) = (exp 



ik- (x(r) - x) + ^ jT dr'^ (x(r')) 



(28) 



The last term in the argument of the exponential accounts for the compress- 
ibility effects in the background turbulence produced by the polarization 
drift. This term has zero average, but as shown below, its correlation with 
the displacements 



U{t) = dr" J dr 1 ( Vi (x(r")) <9 T ,A0 (x(r'))> . (29) 

is not zero and it can play a role in the evolution of the turbulence. 
The ion density is 

Sn'fr, v, t) = -n (x)^ (kyV* ~ up 2 s k 2 ± ) Tf (30) 
and the dispersion relation for test modes on turbulent plasma is 

- fe K„-^)ff=l+Vf^^ (3D 

since the electron density perturbation is the same as in quiescent plasma. 

Thus, the effects of the background turbulence appear in the function of 
time M(r;t) defined in Eq. (|28p . which determines the ion propagator. In 
the case of quiescent plasmas M — 1. This function imbeds all the effects 
of the background turbulence. It depends implicitly on the background po- 
tential through the statistical characteristics of the trajectories (1261) . which 
determine the average, and explicitly through the compressibility term Lj. 

This function and its evolution is estimated in the next sections. 

The polarization drift has an essential role. It destabilizes the drift waves 
but it also has a more complex effect through the background turbulence 
which is a weakly incompressible environment for the test modes due to 
Up. As shown below, the latter effect is important in the strongly nonlinear 
regime. 



3 Test particles in turbulent plasma 

The E x B drift in turbulent plasmas can determine trajectory trapping or 
eddying. It is permanent in the case of static electric fields and appears due 
the Hamiltonian character of the motion with the potential <ft as conserved 
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Hamiltonian. When this motion is weakly perturbed by slow time variation 
of the potential or by other components of the motion, trapping persists for 
finite time intervals determined by the strength of the perturbation |15]-|17j. 

Semi-analytical statistical methods (the decorrelation trajectory method 
DTM [7] and the nested subensemble approach NSA [8]) have been developed 
for the study of test particle stochastic advection. A series of studies ([8] and 
the references there in) have shown that trapping strongly influences the 
statistics of trajectories leading to memory effects, quasi- coherent behavior 
and non-Gaussian distribution. The trapped trajectories have quasi-coherent 
behavior and they form structures similar to fluid vortices. The diffusion 
coefficients decrease due to trapping and their scaling in the parameters of 
the stochastic field is modified [IE]-[22]- Anomalous diffusion regimes and 
even sub-diffusion or super- diffusion can appear due to trajectory trapping. 

DTM and NSA reduce the problem of determining the statistical behavior 
of the stochastic trajectories to the calculation of weighted averages of some 
smooth, deterministic trajectories determined from the Eulerian correlation 
(EC) of the stochastic field. These methods are in agreement with the statis- 
tical consequences of the invariance of the potential. NSA is the development 
of DTM as a systematic expansion that validates DTM and obtains much 
more statistical information. 

The main idea is to study the stochastic equation fl26l) in subensembles 
of realizations of the stochastic field. First the whole set of realizations R 
is separated into subensembles (5*1) that contain all realizations with given 
values of the potential and of the velocity in the starting point of the trajec- 
tories. Then each subensemble (Si) is separated into subensembles (S2) that 
correspond to fixed values of the second derivatives of the potential. By con- 
tinuing this procedure up to an order n, a system of nested subensembles is 
constructed. The stochastic (Eulerian) potential and velocity in a subensem- 
ble are Gaussian fields but nonstationary and nonhomogeneous, with space- 
and time- dependent averages and correlations. The correlations are zero in 
x = 0, t = and increase with distance and time. The average potential 
and the average velocity in a subensemble are determined by the Eulerian 
correlation of the potential by conditional averages (see [8] for details). The 
stochastic equation fT2oT) is studied in each highest-order subensemble. Ne- 
glecting trajectory fluctuations, the average trajectory X(i; Sn) is obtained 
from an equation that has the structure of the equation of motion (126|) . but 
with the stochastic potential replaced by the subensemble average poten- 
tial which is determined by the EC of the potential. This approximation 
is acceptable because it is performed in the subensemble where trajectories 
are similar as they are super-determined. In addition to the necessary and 
sufficient initial condition x(t) = x, supplementary initial conditions are im- 
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posed by the subensemble definition. The strongest condition is the initial 
potential, which is a conserved quantity in the static case and determines 
comparable trajectory sizes in a subensemble. Moreover, the amplitude of 
the velocity fluctuations in subensembles, the source of the trajectory fluctu- 
ations, is zero in the starting point of the trajectories and reaches the value 
corresponding to the whole set of realizations only asymptotically. 

The statistical quantities are obtained as weighted averages of these tra- 
jectories X(t; Sn) or of the functions of these trajectories. The weighting 
factor is the probability that a realization belongs to the subensemble (Sn), 
and is analytically determined. 

NSA is quickly convergent because the mixing of periodic trajectories, 
which characterizes this nonlinear stochastic process, is directly described 
at each order. Results obtained in the first order (DTM) for the diffusion 
coefficient D(t) are essentially not modified in the second order [8]. Second- 
order NSA is important because it provides detailed statistical information on 
trajectories, which contributes to the understanding of the trapping process. 

The EC of the potential is modelled according to the frequency and the 
growth rates of drift modes in quiescent plasma (|T5|) -( TT6|) . This equations 
show that the maximum of 7 corresponds to k±p s = 1, with larger k y (of the 
order of p~ l ) and smaller k x . The growth rate is zero for k y = 0. Thus, the 
spectrum of the drift type turbulence is zero along the line k y = 0, which 
means that the integral of the EC along y must be zero. It has two peaks 
at k y = ±k Q , k x = 0. A simple correlation with these properties that also 
includes anisotropy is 



where (3 is the amplitude of the fluctuations of the potential, 6 = X^/Xy, the 
distances are normalized with A = X y and k with A -1 . This correlation also 
accounts for the existing of dominant waves in the stochastic potential. This 
correlation has negative values along y (one negative minimum for ko = 0, 
and oscillatory behavior for large ko). 

Starting from the EC (!32|) . the statistics of ion trajectories is determined, 
then the averaged ion propagator, the growth rate and the frequency of the 
drift modes are estimated. These quantities show how the correlation evolves. 

The distribution of displacements obtained from Eqs. f l26]) strongly de- 
pends on the ordering of the characteristic tines of the stochastic process. 




(32) 



f(y) 



cos k y - y 



sin k y 
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The longest characteristic time is the parallel time of the ions t\u = Xn/v^i, 
which actually leads to negligible parallel motion. The motion of the po- 
tential along y with the diamagnetic velocity defines the diamagnetic time 
t* = X y /V* e . The correlation time of the potential r c is the characteristic 
time for the change of the shape of the potential r c = 7 -1 , which is larger 
than the diamagnetic time (r c > r*). The time of flight (or eddying time) is 
defined by tji = X x /V x = BX x X y /(3. Trajectory trapping appears when tji is 
the smallest of these characteristic times. The trapping parameter for drift 
turbulence is 

K* = n/r fl = V y /V* e = 0/(X x V* e B) 

and this process appears when K* > 1. 

The average propagator fT271) is calculated in the next sections using sim- 
plified models that include the main characteristics of the probability of dis- 
placements. It is so possible to capture the complicated nonlinear effects 
of strong turbulence in rather simple analytical expressions. The results 
can easily be improved by taking into account the statistics of test particles 
obtained with the nested subensemble method. This significantly more com- 
plicated approach that relies on numerical calculation of the averages does 
not change qualitatively the results. We consider that the simple analytical 
expressions derived in the next sections give a more clear image on the com- 
plex nonlinear processes that appear in the drift turbulence evolution beyond 
the quasi-linear stage. 

4 Ion diffusion and damping of large k modes 

For small amplitude of the background potential, the time of flight is larger 
than the decorrelation time and trapping does not appear. The smallest 
characteristic time that influence the ion motion corresponds to the poten- 
tial motion with the diamagnetic velocity, which travels over the correlation 
length faster than the ions. The characteristic times ordering is 

T|| e < T* <T C < Tfi < T\\i. (33) 

This corresponds to small amplitude background turbulence with (3/BX x <C 
V* e , which does not produce trajectory trapping. The EC of the turbulence 
is given by (I32p that corresponds to the growth rates of the drift modes. 

The diffusion coefficients determined for this quasilinear regime for a 
frozen potential that moves with V* e and has the EC (1321) are completely 
different from those obtained for an EC that is a decaying function of the 



13 



distance. In the latter case, the transport along the average velocity is dif- 
fusive with D y = VyT* = VyXy/V^ and the perpendicular transport is sub- 
diffusive with the decay of the diffusion coefficient given by the EC of the 
potential D x (t) ~ \d y E(V* e t)\ /V* e .The EC (1321) determines similar subdif- 
fusion along x and subdiffusion along y also. The potential change on the 
longer time scale r c leads to diffusive transport but with very small diffusion 
coefficients on both directions, much smaller than for a monotonically decay- 
ing potential. The dependence of the diffusion coefficients on r c is different 
on the two directions: D x increases while D y decreases when r c increases. 
The distribution of the displacements is Gaussian, as is the velocity. 

The function fl28l) calculated with the Gaussian distribution is given by 



M(t; t) = exp 



kj (fo(r) - x,) 2 ) ' 
+ ikiLi 



(34) 



where the quadratic term in the divergence of the polarization drift was 
neglected. The quantity Liij) defined in Eq. (1291) is a Lagrangian correlation 
that has the dimension of a length and represents the correlation of the 
compressibility term with ion trajectories in the background turbulence. 

The trajectories (the characteristics) have displacements until decorrela- 
tion by potential motions that are small compared to the correlation length 
and can be neglected. This strongly simplifies the estimation of the com- 
pressibility term Lj because the Lagrangian correlation in (129]) can be ap- 
proximated by the corresponding Eulerian correlations. The correlation in 
L x is 

K(x(r")) ^A0(x(r'))} = (d x <j> (x - V* e (t - r")) d T ,A0 (x - V* e (t - r'))) 

= -Ld T „d Tl AE(V* e (T"-T')) 

' *e 

which integrated over time gives 

L,(r) = (AE (0) - AE (V« \t - r|)) , (35) 

D V* e 

which is finite because the Laplacean of the EC is symmetrical and finite in 
zero. This function of time starts from zero at time r = t and saturates in a 
time t— r of the order of r* at the value 2AE (0) / BV* e . The other component 
L y is zero because it contains one or three derivatives at x for x = 0, which 
are zero. 



14 



Thus in the quasilinear case 



M = exp 
= exp 



-kfDiit - r) + %k x 



-kfDi(t - r) - i2ki 



eB 3 V M 
V 2 



(36) 



which shows that the compressibility determines a constant term propor- 
tional with the square of the amplitude of the stochastic velocity determined 
by the background turbulence. This is very small due to and to the 
condition V <C V* e imposed for this regime. 
The propagator is 



IT = — exp 
and the ion density perturbation 



-i2fci 



V 2 



(37) 



^(x,,,t)=n (a;)^ A; ^ e -^ A;i 



where 



exp 



V 



2 1 



to + ikfDi 
V 2 

1 - i1k x 



(3? 



(39) 



This factor produced by the polarization drift is very close to unity. 
The dispersion relations is 



fc y K e - up 2 s k\ _ Fk 
co + rk 2 D t A - l + l \j2 



co kyV^ e 



\k z \v Te 



and its solution is 



A', 



CO 



1 + k 



7 = 7 o; (k y - co) — k { Di — 2k\co 



PsC s 



(40) 



(41) 



where Di = Di/(p s V* e ). Thus the frequency is not influenced by the back- 
ground turbulence when it has small amplitude. The latter determines new 
terms in the growth rate. 

The main effect is produced by ion trajectory diffusion in the background 
turbulence, which determine the damping of the modes with large k (the sec- 
ond term in 7). This is the well known result of Dupree, which is recovered 
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here. An additional effect is obtained from the compressibility, which deter- 
mines the asymmetrical increase of the modes (contribution to the growth 
of the modes with kik2 < and to the decay of those with kik2 > 0). This 
effect is small because the ratio of the third term to the first term in Eq. ( 140]) 
is of the order (l/^r|| e ) (V 2 /V? e ) < 1. 

The frequency and the growth rate of the modes on turbulent plasma 
with small amplitude show that the most unstable modes are, as in quiescent 
plasma, those with u = k y V* e /2, which have wave numbers of the order of 
Pj 1 — !)• Their interaction with the background turbulence is weak, 
except for the large k modes, which decay and eventually are damped due to 
ion diffusion. 

This means that EC of the turbulence evolves in this quasilinear stage by 
increasing its amplitude (3 without major change of its shape. 

5 Trajectory structures and large scale cor- 
relations 

The increase of the amplitude (3 of the stochastic potential determines the 
change of the ordering of the characteristic times (133|) . When V becomes 
larger than the velocity of the potential V* e (or (3 > V* e \ y B), the time of 
flight is smaller that the decorrelation time r* 

T|| < Tfi < T* < r c < rj. (42) 

In these conditions ion trapping or eddying appears. As we have shown 
this strongly influences the statistics of trajectories. The distribution of the 
trajectories is not more Gaussian due to trapped trajectories that form quasi- 
coherent structures. At this stage the trapping is weak in the sense that the 
fraction of trapped trajectories n tr is much smaller than the fraction rif of 
free trajectories. 

The probability of displacements P(x, t) was obtained using the nested 
subensemble method. It has a pronounced peaked shape compared to the 
Gaussian probability. The contribution of the trapped particles to the prob- 
ability P(x, t) at time t appears in the narrow peak in r = 0, which remains 
invariant after formation. The time invariance of this central part of the 
probability is due to the vortical motion of the trapped particles. It shows 
that vortical structures appear due to trapping. The size of the vortical struc- 
ture depends on the strength of the decorrelation mechanism that releases 
the trapped particles. The free particles, for which the motion is essentially 
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radial, determines the larger distance part of the probability, which continu- 
ously extends. The probability of displacements until decorrelation (r < r*) 
is modeled by 

P(x, y, r) = n tr G(x; S) + n/G(x; S') (43) 

where G(x; S) is the 2-dimensional Gaussian distribution with dispersion 
S = (S x ,S y ). The first term describes the trapped trajectories. We have 
considered for simplicity their distribution as a Gaussian function but with 
small (fixed) dispersion that represents the average size of the structures. 
The shape of this function does not change much the estimations. The free 
trajectories are described by the second term in Eq. (j4"3"j) . They have dis- 
persion that grows linearly in time: (r) = Si + 2Di (t — r) , i — x,y. The 
initial value S' { (t) = Si is an effect of trapping. It essentially means that the 
trajectories are spread over a surface of the order of the size of the trajectory 
structures when they are released by a decorrelation mechanism. 

The compressibility term L, will be calculated in the next section. As 
shown there, it is negligible in these conditions. 

The average propagator is 

= Wft^ ( 44 ) 

where the factor T is determined by the average size y/S~i of the trapped 
trajectory structures 

^ = e*pf-i*?S<Y (45) 
The dispersion relation fl3~Tl) is 

(k y V* e - up 2 s k\) . Fku - k y V* e 



- F=1 + H r_ z _ J r I e (4g) 



u + ikfDi ' V 2 \k z \v T , 

with the solution _ 

57 = (47) 
1 + Tk^ 

7 = 7 cU (k y - uJ) - fc- A. (48) 

Thus the effect of trajectory trapping appears in the frequency while the 
growth rate is not modified (the small compressibility term was neglected). 
The trajectory trapping determines the decrease of the frequency. The value 
of k± corresponding to the maximum of 7 (at u7 = k y /2) is displaced from 
values k± = 1 to k± = 1/s where s = S x + S y /p s . For s ^> 1 (Si ^ p 2 s ) , 
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Figure 1: The frequencies (dashed curves) and the growth rates of drift modes 
for several values of the size of trajectory structures 

the frequency (1471) is UJ = J-'ky, which shows that the finite Larmor radius 
effects become negligible in the presence of trapping when the size of the 
vortical structures is larger than p s . 

The growth rate of the drift modes for different values of the size s of 
the trajectory structures is shown in Figure 1. This shows that the range of 
the wave numbers for the unstable test modes is displaced toward small kj_. 
This determines a further increase of the size of trajectory structures. The 
diffusive term in (T4"8"l) lead to the damping of the modes with k±_ 3> 1/s. 

This means that EC of the turbulence evolves in this weakly nonlinear 
stage by increasing its amplitude (3 and by increasing the correlation length 
and the average wave number k (1321) . The vortical structures of ion trajec- 
tories produced by trapping determine this processes. 

6 Ion flows and turbulence attenuation 

The evolution of the potential determines the increase of the fraction of 
trapped ions. This determines an average flux of the trapped particles that 
move with the potential n tr V* e . As the E x B drift has zero divergence, the 
probability of the Lagrangian velocity is time invariant, i. e. it is the same 
with the probability of the Eulerian velocity. The average Eulerian velocity 
is zero and thus the flux of the trapped ions that move with the potential 
has to be compensated by a flux of the free particles. These particles have 
an average motion in the opposite direction with a velocity Vf such that 

n tr V* e + n f V f = 0. (49) 
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The NSA shows that the probability of the displacements splits in two com- 
ponents that move in opposite direction. The peak of trapped ions moves 
with the velocity V* e while the free ions move in the opposite direction with 
a velocity that increases with the increase of the fraction of trapped ions. 
The distribution of trapped ions is almost frozen while that of free ions has 
a continuously growing width. Thus, opposite ion flows are generated by the 
moving potential in the presence of trapping. A simple approximation that 
includes the main features of the distribution is 

P(x,y,t) = n tr G(x,y- V* e t; S x , S y ) +n f G(x,y- V f t; S' x , S'). (50) 



where two Gaussian functions were used as in Eq. (1431) but taking into 
account the ion flows fT49]) . The average velocity of the free ions is Vf = 
—nV* e , where n = n tr /nf. The separation of the distribution and the exis- 
tence of ion flows in drift type turbulence are confirmed by numerical simu- 
lations [14J. 

The compressibility term Lj (T2"9"l) is determined for trapped and free par- 
ticles with the DTM. It can be written as 

L x = - J^j jf dr" jf dr' {d y cj> (x(r")) <9 y A0 (x(r'))> . 

A nonzero value is obtained for L x essentially because the potential and the 
Laplacean of the potential are correlated in a stochastic field: the maxima 
of the potential statistically coincide with the minima of the Laplacean. It 
follows that the derivatives of these functions along the same directions are 
also correlated and the derivatives on different directions are not. Indeed, 
one obtains from Eq. 



(cW>(x) <9 y A0(x)} = -d x d v AE(0) = (51) 

?(ri\ = - 

A 4 



(<9 y 0(x) d y A<f>(x)) = -d 2 y AE(0) = (52) 



where 



k 2 n + 3 



+ 15 + 10 kl + k%) (53) 



The corresponding Lagrangian correlations are determined with the DTM. 
The component L y is zero for both trapped and free ions due to the symme- 
try of the correlation (1511) . L x is different for free and trapped trajectories. 
In the first case there is a relative motion of the free ions and of the poten- 
tial with the velocity V* e /rif (obtained from (j49| ). which determines rapid 
decorrelation. The estimation done for the quasilinear turbulence also holds 
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in these conditions and L x for free ions is approximated by Eq. (1551) with V* e 
replaced by V* e /rif, and thus it is negligible. For the trapped trajectories, 
L x is much larger than in the quasilinear case. It can be determined with 
the DTM, which shows that the following estimation holds. 



L x {r, t) 



a 5 2 



dr" / dr'C x {\r" -r'\) 



C x (\r" - r'\) = (d y ^ (x(r") - V* e (t - r")) d y A0 (x(r') - V-(t - r'))> . 



The integral is calculated by changing the integration variable r" to 9 

2V*, 



T " -T 



t-T 



which integrated by parts gives 

L x (r,t) = - 



d6C x (\6\)(t 



2Ke 

as 2 



t-T 



d9 J{6) 



where 



drC r 



T 



This function saturates after the decorrelation time at J (9) = C x (0) Ta- 
in these conditions of strong turbulence, the decorrelation time for the C x is, 
as for the correlation of the stochastic E x B velocity, the time of flight r/j. 
Thus 



L x {t, t) = a(t-r), a = 2<9 2 A£(0) 



TflV M 



(54) 



where a has the dimension of a velocity. This type of Lagrangian correlation 
corresponds to the decorrelation by mixing. 

Thus, the compressibility of the background turbulence is correlated with 
the displacements along x and influences the propagator for the trapped 
ions. The function M ( 128|) that accounts for the effects of the background 
turbulence is determined using Eqs. (150]) and ([51 



M 



nf exp 
n tr exp 



- k l{S^ + 2D l {t-r))- l k y V f {t-r) 



+ 



(55) 
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The average propagator (1271) becomes in these conditions 



n = —T 



+ 



n tr 



u — kyVf u — kyV^ + k x a 



- ikfDi 



n f 



{u - kyVf) 

where the fact that the diffusion term is small was used. 
The dispersion relation is 



(56) 



T (kyV* e - uk 2 p 2 s ) 



n 



f 



n tr 



k 2 Dm 



Ul 



kyV^ e + k x a (u - k„V, 



y v f) 



l + i 



7T U 



h V 



2 \k z \v T e 



(57) 



The compressibility term k x a competes with the term u — k y V* e that is 
of the order k y V* e . Its effect is measured by the ratio 



k x a 
h V 



2c 



a 



~Kp s v 



(5? 



where c is the normalized correlation (|5U|) . For a typical plasma with c 

-3 



10 6 m/s, 13 = 150V, B = 3T, A = 10~ 2 m, p s = 2 . ,„„ , , , s - 

5 10 -3 . The correlation c > 20, which leads to r = 0.05 for k y = k x . It means 
that the compressibility term can be neglected for modes with k y larger or 
of the order of k x and that it is important for modes with k y — > 0. 

The solution of the dispersion relation is determined for the two cases: 
k y > k x that corresponds to drift modes and k y = 0, which are modes with 
different characteristics called zonal flow modes. 



6.1 Drift modes 

Neglecting the compressibility term (a = 0), the real part of the dispersion 
relation ([57]) is 



T [k y - Udk]^ [cu d + (n-l) k y ] = (u d + nk y ) (u d - k y ) , (59) 

where the suscript was introduced in u d to account for drift type modes. 
Taking u d = U + 5(n, k), where ZJ = Tk y j ^1 + Fk 2 ^j is the frequency (HTj) 
in the absence of ion flows (n = 0), the equation becomes 

— i — — -i Tik„. 

5 p + aJo + (n - 1) ky] = y —, (60) 

1 + Fk ± 
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where Doppler shift with the velocity of the trapped and free ions appears in 
the parenthesis in the left side term, uj —k y Vf—k y V* e = k y V* e \uJq + (n — 1) k y \ 
The solutions are the two intersections of the parabola in the left hand side 

of Eq. f )60|) with the horizontal line at the level nk y / ^1 + Fk 2 ^ . The con- 
dition 5 = for n = chooses one of the two solutions, namely on the 
branch of the parabola that passes through the origin. This leads to 5 < if 
uJq + (n — 1) k y < and 5 > if cJ + (n — 1) k y > 0. The absolute value of 
5 is an increasing function of n in both cases, and thus, as n increases, the 
decrease of uJd is produced by the ion flows in the first case and the increase 
of ujd in the second case. The solution of the dispersion relation fl59l) is 



Ud= 2 



(jj + (1 — n) k y + sg 



\ 



in 



l)ky] 



2 4:nkl 
+ — ' 



l + Tk] 



where sg = sign [wo + (n — 1) k y ~j . A jump appears in the frequency at the 
value of kd determined by the equation uJ = (I — n) k y , which has solution if 
n < 1. kd is an increasing function of n.In the limit of large k, uJd — > —nk y if 
n < 1 and ZJd — > k y if n > 1 . This means that the asymptotic phase velocity is 
in the ion diamagnetic direction when n < 1 and in the electron diamagnetic 
direction when trapping is stronger such that n > 1. 

The imaginary part (obtained using the condition 7 d uJd) is 



2 \ 2 

1 - uj d k ± ) f+ + k ± 



■=- + 



n tr 



uj d + nky uj d - k 



_ _ 2 _ F f 1 - v d k ± ) 

7o (k y - uj d ) - n f k i Di— — / 

(ujd + nky) 



where 



(uJ d + nky) (uJ d - k y ) 
Using the solution can be written as 



+ 



n tr 



1- \2 



> 0. 



7o (uj d + nky) [(1 -n)k y - uj d ] - n/A; f A _ , : 

7d = — ~ — — ~=^\ =2— l k v - ^d) 



[(1 — n)k y — uj d ] ( 1 + Fk ± J + nk 



-2 



(61) 



The growth rates are presented in Figure 2a. The jump in the frequency 
corresponds to a jump in the growth rate from negative (at small k) to 
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positive values. Thus the large wave-length modes are stabilized by the 
ion flows, while the small wave-length modes still grow. This leads to the 
increase of the amplitude of turbulence accompanied by the decrease of its 
correlation length. Both effects contribute to the increase of the ratio r*/ r/; = 
P/(V* e \ y B), which leads to the increase of n. Consequently, the increase of 
the amplitude of the turbulence continues as well as the damping of small 
k modes and the decrease of the correlation length. The process is stopped 
when n > 1 because the growth rate becomes negative for the whole wave 
number range. This determines the decrease of n and the attenuation of the 
ion flows. 




Figure 2: Growth rate dependence on the fraction of trapped ions n and on 
the size S of trajectory structures for drift modes (left panel) and zonal flow 
modes (right panel). The curves with the same label correspond to the same 
parameters. 



6.2 Zonal flow modes (k y = 0) 

The compressibility term becomes dominant for k y = 0. As shown below it 
determines unstable modes completely different of the drift modes: with k y = 
0, k x 7^ and very small frequency u z f, much smaller than the diamagnetic 
frequency. These are zonal flow modes. 



The dispersion relation for k v = is 



-UzfTklpl 



tUzf + ik^Di u z f + k x a 



2 k?_ 



VTe 



For 7 2 j, kfDi <ti u z f, the real part is 



-u zf J r k x p s (u zf + n f k x a) =u zf (u zf + k x a) 
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with the solution 

Uzf = 

where 

_ a 

The imaginary part is 

pK 4 2 71/ 

= \ -T\\ e u zf -k x D x F Ps 

V 2 w 2 / 

with the solution 

TltrJ~ k x 

(63) 

These unstable zonal flow modes are consequence of trapping combined with 
the polarization drift. One can see that when trapping is negligible (n tr — 0), 
u z f = —k x a and j z f = 0, and that u z j, j z f = for a = 0. As seen in Figure 
2b, the growth rate of zonal flows is positive and it increases with the increase 
of n tr and with the decrease of the size of the trajectory structures. 

Thus we have shown that zonal flow modes are generated due to ion flows. 
They have positive growth rates and small frequencies (typically ten times 
smaller than the diamagnetic frequency). 

In conclusion, the ion flows produced by ion trapping in the moving po- 
tential determine two parallel effects: nonlinear damping of drift modes and 
generation of zonal flow modes. There is no causality relation between zonal 
flow modes and drift turbulence damping. Both effects are generated self- 
consist ently in the nonlinear evolution of drift turbulence. The influence 
produced by the zonal flows on the drift type modes is only indirect, through 
the diffusive damping. Zonal flows modify the configuration of the potential 
and consequently the correlation of the turbulence by introducing compo- 
nents with k y = in the spectrum, which change the shape of the EC. The 
integral of the EC along y becomes finite. This determines a rather strong 
increase of the diffusion coefficient D y , which contribute to the decay of the 
drift type turbulence. Zonal flow modes also modify the diffusion along x 



1 + 



(62) 



■ = 2c 



a 



~p s V 



kl + 3 A c s 




n tr 



n tr 



(u zf + k x a) 2 ) u z f u zf + k x a 



-r \ 2 



7o (k x a) 



-2\' 

(l + ^) 2 



- n s 7k x D x 
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but in the opposite sense: a strong decrease appears due to the decrease of 
X x produced by the zonal flow component of the turbulence. The drift tur- 
bulence decay determines the decrease of n tr , which reduces the growth rate 
of zonal flow modes. 

There is time correlation between the maximum growth rate of zonal 
flow modes and the damping of the drift modes, as observed in experiments 
and in the numerical simulations. This time correlation can be deduced 
from Figures 2 (a, b), which show the growth rates of the drift and zonal 
flow modes for the same set of parameters n, S. When the drift turbulence 
amplitude increases determining the increase of the fraction n of trapped 
ions (curves 1-5) the growth rate of the zonal flows is very small. When 
the growth rate of the drift modes becomes negative and the size of the ion 
trajectory structures decreases (curves 6, 7 in Figure 2) the growth rate of the 
zonal flow modes strongly increases and, at small S, it becomes comparable 
with that of the drift modes. The attenuation of the drift turbulence that 
leads to the decrease of n determines the decrease of the growth rate of the 
zonal flow modes (curves 8, 9 in Figure 2 b). 

We note that the effect of ion diffusion on the zonal flow modes and of 
drift modes is different. Ion diffusion determines the damping of the drift 
modes with large k. This effect is well known and appears on the whole 
evolution of the drift modes, from the quasilinear to the strongly nonlinear 
stage characterized by the existence of ion flows. In the case of zonal flow 
modes, Eq. (??) shows that the effect of ion diffusion on large k x modes is 
negligible due to the factor J 7 that goes to zero. The diffusion has strongest 
effect on the modes with k x = y/l/S x , and it can even determine negative 7 
at small k x . 

The evolution of the turbulence in this strongly nonlinear stage is rather 
complex. The ion flows determine first the increase of the amplitude of 
the turbulence and the decrease of the correlation length X y , then a major 
modification of the EC shape appears due to generation of zonal flow modes 
and eventually the turbulence is strongly attenuated. The new component 
that accounts for the zonal flow modes adds to the EC a function with a slow 
decay along y (k y = 0), with zero integral along x (because 7 = for k x = 0) 
and with amplitude growing with the increase of n. It essentially determines 
in the total EC the increase of \ y and of the amplitude (3 and the decrease 
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7 Summary and conclusions 



Test modes on turbulent plasmas were studied taking into account the pro- 
cess of ion trajectory trapping in the structure of the background potential. 
The case of drift turbulence was considered, and the frequency u> and the 
growth rate 7 were determined as functions of the statistical properties of 
the background turbulence. The main characteristics of the evolution of the 
turbulence were deduced from 7 and u. 

A different physical perspective on the nonlinear evolution of drift turbu- 
lence is obtained. 

The main role is played by the trapping of the ions in the stochastic po- 
tential that moves with the diamagnetic velocity. The fraction of trapped 
ions is a function of = (3/ (\ x V* e B). Trapping appears for > 1 and the 
fraction of trapped trajectories rit r increases with the increase of this param- 
eter. Trapped trajectories form quasi- coherent structures, which determine 
non-Gaussian distribution of the displacements. At larger amplitudes, the 
moving potential determines ion local flows: the trapped ions move with the 
potential while the other ions drift in the opposite direction. These opposite 
(zonal) flows compensate such that the average ion velocity is zero, but they 
determine the spliting of the distribution of displacements. 

The first effect of trapping in turbulence evolution appears when the 
fraction of trapped ions is small and the ion flows are negligible. The trapped 
ions determine the evolution of the turbulence toward large wave lengths (the 
inverse cascade). 

The influence of the ion flows produced by the moving potential appears 
later in the evolution of the turbulence when the fraction of trapped ions 
is comparable whith that of free ions. The ion flows determine first the 
damping of the small k±_ drift modes and eventually the damping of the 
drift modes with any k±. In the same time, trapping in connection with 
the compressibility produced by the polarization drift in the background 
turbulence determine transitory zonal flow modes (with k y = and uj <C 
kyV* e ). Thus, in this perspective, there is no causality connection between 
the damping of the drift turbulence and the zonal flow modes. Both processes 
are produced by ion trapping in the moving potential. 

During the decay of the turbulence trapping is maintained because the de- 
crease of the amplitude is compensated by decrease of the correlation length, 
leading to small variations of the parameter K*. The damping is first pro- 
duced for the small wave numbers and thus there is a transitory domination 
by large wave numbers in the spectrum. This determines a hysteresis process: 
the growth and the decay of the turbulence is produced on different paths. 
Large scales are generated at the development (increase) of the turbulence 
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and the correlation length increases. This determines only a slow increase of 
the parameter and thus a slow increase of the number of trapped trajec- 
tories n^.Then, the ion flows become important and produce further increase 
of the larger part of the wave number spectrum accompanied by the decay 
of the small k modes. Thus, the amplitude of the potential increases and 
the correlation length decrease. This determines a much faster increase of 
and of the relative number of trapped ions. When the latter is n tr = 0.5, 
the growth rates of all drift modes are negative and the amplitude of the 
turbulence decays leading to the decrease of rit r . This process continues until 
the ion flows become negligible {n tr A closed evolution curve in the 

(j3, A) space is described by the turbulence, which remains in the nonlin- 
ear stage characterized by trapping and oscillates between weak and strong 
trapping. 

The zonal flows appear due to ion flows as well as turbulence nonlinear 
damping. The damping process is not determined by the zonal flows. There 
is however an influence produced by the zonal flows on the drift type modes. 
It is not direct but through the diffusive damping. Zonal flows modify the 
correlation of the turbulence by introducing components with k y = in the 
spectrum. This changes the shape of the EC, which has no more zero integral 
along the y direction. This determines the increase of the diffusion D y , which 
contribute to the decay drift type turbulence but not of the zonal flow modes. 

The diffusion coefficient along x axis also oscillates during the evolution 
of the turbulence. Its maximum is correlated with the period when the 
correlation lengths along x is large, which appears before the development of 
the ion flows. 

The zonal flow modes increase and drift turbulence decay appear to be 
temporally correlated, as observed in experiments and simulations. The char- 
acteristic time At for turbulence and transport oscillations can be estimated 
as the inverse of the growth rates, which are of the order of 5 10~ 3 c s /L„. 
One obtains for L„/i? = 5, At = 40 Ro/c s , in agreement with the numerical 
simulations. 

This work was supported by CNCSIS-UEFISCDI, project number PNII 
- IDEI 1104/2008. 
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